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We show how many-body ground state entanglement information may be extracted from sub- 
system energy measurements at zero temperature. Generically, the larger the measured energy 
fluctuations are, the larger the entanglement is. Examples are given with the two-state system 
and the harmonic oscillator. Comparisons made with recent qubit experiments show this type of 
measurement provides another method to quantify entanglement with the environment. 
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A many-body quantum system is cooled to zero tem- 
perature so that it is forced into its overall nondegen- 
erate ground state. We discuss the measurement of a 
sub-system Hamiltonian and demonstrate that it can be 
found in an excited state with a probability that depends 
on the coupling to its environment. This non- intuitive re- 
sult is a pure quantum phenomenon: it is a consequence 
of entanglement 0] of the sub-system with the environ- 
ment. In fact, we demonstrate that knowledge of the 
probability to find the system in an excited state can be 
used to determine the degree of entanglement of the sub- 
system and bath. Consequently, simple systems with well 
known isolated quantum mechanical properties (such as 
the two-state system and harmonic oscillator) become 
"entanglement-meters" . 

There is growing interest of ground-state entangle- 
ment in condensed matter physics. Theoretical works on 
ground state entanglement have addressed entropy scal- 
ing in harmonic networks 0, spin-spin entanglement in 
uantum spin chains 0] and quantum phase transitions 
""J. Entanglement properties of the ground state are 
also essential in the field of adiabatic quantum comput- 
ing 6] . Furthermore it is interesting to link other ground 
state properties like the persistent current of small meso- 
scopic rings Li l | g|h or of doubly connected Cooper pah- 
boxes 0, fiol lllj. or the occupation of resonant states 
[l2| to the zero-temperature entanglement energetics. 

We consider a general Hamiltonian H = H s +H c +He, 
that couples (c) the system (s) we are interested in to 
a quantum environment (E) such as a network of har- 
monic oscillators 0] . The lowest energy separable state 
is |fi) = |0) s |0)_e, where |0){ Sj £} are the lowest uncoupled 
energy state of both systems. However, if the system 
Hamiltonian and the total Hamiltonian do not commute 
(which is the generic situation), then |fi) is not an energy 
eigenstate of the total Hamiltonian. Thus, there must be 
a lower energy eigenstate (|0)) of the total Hamiltonian 
which is by definition an entangled state. Because time 
evolution is governed by the full Hamiltonian, the ground 
state expectation of any operator with no explicit time 
dependence will have no time evolution, insuring that 
any measurement is static in time. This situation is in 
contrast to the usual starting point of assuming that the 



initial state is a separable state and studying how it be- 
comes entangled. The reduced density operator of the 
system is given by tracing out the environmental degrees 
of freedom, p = Tr^|0)(0|. Assuming the full state of 
the whole system is pure, the reduced density matrix 
contains all accessible system information, including en- 
tanglement of the system with its environment. Because 
repeated measurements of H s will give different energies 
as the sub-system is not in an energy eigenstate, we are 
interested in a complete description of the statistical en- 
ergy fluctuations. These fluctuations may be described in 
two equivalent ways. The first way is to find the diagonal 
density matrix elements in the basis where H s is diago- 
nal. These elements represent the probability to mea- 
sure a particular excited state of H s . A second way is to 
find all energy cumulants. A cumulant of arbitrary order 
may be calculated from the sub-system energy generating 
function, Z(\) — (exp(— x-ff s )) (as always, (O) = TrpO) 
so that the n th energy cumulant is given by 



(i) 



These cumulants give information about the measured 
energy distribution around the average. 

Before proceeding to calculate these energy fluctua- 
tions, we ask a general question about entanglement. 
Given the energy distribution function (the diagonal ma- 
trix elements of the density matrix only), can anything be 
said in general about the purity or entropy of the state? 
Surprisingly, because we are given the additional infor- 
mation that we are at zero temperature, the answer is 
yes. If we ever measure the sub-system's energy and find 
an excited energy, then we know the state is entangled. 
Though this statement alone links energy fluctuations 
with entanglement, a further quantitative statement may 
be made in the weak coupling limit. The reason for this 
is the following: the assumptions exponentially suppress 
higher states, so to first order in the coupling constant, we 
can consider a two-state system where the density matrix 

has the form p = 1 — ap, p++ = ap, p^ = p*_, = ac. 

For vanishing coupling constant a = 0, this just gives 
the density matrix for the separable state. The linear 
dependence of p on a holds to first order for the model 



systems considered below and is the entanglement contri- 
bution. If one measures the diagonal elements of p, one 
obtains and p ++ as the probability to be measured 
in the ground or excited state (because a is small, there 
is only a small probability to find the sub-system in the 
upper state). If we now diagonalize p, the eigenvalues 
are \± = {1 — pa, pa} + 0(a 2 ). To first order in a, the 
eigenvalues arc the diagonal matrix elements, so wc may 
(to a good approximation) write the purity or entropy in 
terms of these probabilities even if the energy difference 
remains unknown. 

The Qubit. Let us now first evaluate the energy fluctu- 
ations of a qubit, a two-state system. The most general 
(trace 1) spin density matrix is p = (1 + (o~ x )o- x +(a v )a v + 
(o"z)cz)/2. A simple measure of the entanglement is given 
by the purity, Trp 2 = (1/2)(1 + X 2 + Y 2 + Z 2 ), where 
Xi = (<Ji). It is well known that (X,Y,Z) form coor- 
dinates in the Block sphere. Purity lies at the surface 
where X 2 + Y 2 + Z 2 = 1, whereas corruption lies deep 
in the middle. 

We take the system Hamiltonian 0] to be H s = 
(e/2)a z + (A/2)a x . Introducing the frequency ft — 
Ve 2 + A 2 /h and using the identity e ~ l ^ n ' CT = Icos| — 
in ■ a sin 4 with (3 — xfi, it is straightforward to show 



Z(i X ) = cos(ft X /2) 



,sm(0 X /2), 

1 — n — {e{a * 



AK». (2) 



The energy probability distribution may be easily found 
by Fourier transforming Eq. (J2J), or by tracing in the 
diagonal basis of the system Hamiltonian. The an- 
swer may be expressed with only the average energy, 
(H s ) — |(cr 2 ) + y(c x ) as a sum of delta functions at 
the system energies ±M7/2 with weights of the diagonal 
density matrix elements, 



P++,- 



(1/2) 



1 ± 



or. 
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(3) 



Clearly, if the spin is isolated from the environment, 
(H s ) = — fofl/2, the ground state energy, the probability 
weight to be in an excited state vanishes. This distribu- 
tion may also be found from knowledge of the isolated 
eigenenergies, the fact that (H s ) — J2jEjPjj, and that 
Tip = 1. This later argument may be extended to n-state 
systems given the first n — 1 moments of the Hamiltonian 
and the n eigenenergies. 

Connection with Real Qubits. The probability weights 
depend on the energy parameters e and A, and the expec- 
tation values of the Pauli matrices. For real qubits pro- 
duced in the lab, these will depend on the environment 
(lif. Often, we can link the basic phenomena we have 
been describing to physical measurements other than en- 
ergy. For example, in a mesoscopic ring threaded by an 
Aharonov-Bohm flux <3>, or for the Cooper pair box, the 
tunneling matrix element depends on flux A($). The 
operator one is interested in measuring is given by the 
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Coupling constant a 

FIG. 1: Energy probabilities p up and Pdown for the spin-boson 
problem. With increasing coupling to the environment it it 
is more likely to measure the qubit in the excited state. For 
the symmetric case (e = 0), we use the Bethe ansatz solu- 
tion, while for the general case, we use a perturbative solution 
which is only valid for large e or large a. 



normalized persistent current 0, or the expectation of 
a x (for the symmetric case of e = 0), 



mi 



(4) 



where Iq(Q) is the uncoupled value of the persistent 
current. A common model for environmental effects is 
given by coupling the two-state system to a series of har- 
monic oscillators, the spin-boson model 00,0,0. In 
Fig.(l) we have plotted the upper and lower occupation 
probabilities for the spin-boson model as a function of 
the coupling constant a. For the symmetric case (e = 0), 
we have used the Bethe ansatz solution 0, 0, while 
for e finite, we have used the perturbative solution in 
A/u) c which is valid only for larger a or e 7]. Thus the 
plot is cut off at a small a. A computational approach 
calculating the expectation values of the Pauli matrices 
over the whole parameter range was given in Ref. [l5| . 
For e = 0, Eq. J3Jl determines also the expectation value 
of H s , and thus the weights of the energy distribution 
function, Eq. 0, are directly related to the persistent 
current. 

Experiments are always carried out at finite tempera- 
ture, and it is important to demonstrate that there exists 
a cross-over temperature to the quantum behavior dis- 
cussed here. In the low temperature limit, the thermal 



occupation probability is p u% 



-(E 2 - El )/kT_ j n the 



weak coupling limit for the symmetric spin boson prob- 
lem, the probability to measure the excited state scales 
as Pup = — alog(A/w c ) 20|. Setting these factors equal 
and solving for T* yields 



kT* 



E 2 — Ei 
log(alog^) 



(5) 
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Since T* scales as the inverse logarithm of the coupling 
constant, it is experimentally possible to reach a regime 
where thermal excitation is negligible. 

As an order of magnitude estimate, we compare with 
the Cooper pair box ({J [Toj| which is among the most en- 
vironmentally isolated solid state qubits yQJ. From 
which found aQ~ 10 4 , we estimate the quantum prob- 
ability for the box to be measured in the excited state as 
Pup ~ 10~ 3 — 10 -4 , which is of same order or larger than 
the thermal probability, pth ~ 10~ 4 . Experimentally, p up 
and pth may be confused by fitting data with an effec- 
tive temperature, p t h — exp(— f3 e gH s ) |2l|. However, one 
may distinguish true thermal behavior from the effect de- 
scribed here because p up and pth depend differently on 
tunable system parameters such as A. In fact, /3 e ff is an 
entanglement measure. This chain of reasoning may be 
inverted to provide an estimate for Q given only p up . 

The Harmonic Oscillator. We now consider the en- 
tanglement energetics of a harmonic oscillator, H s = 
p 2 /(2m) + (l/2)mu> 2 q 2 . Since there are an infinite num- 
ber of states, the problem is harder. To simplify our task, 
we assume a linear coupling with a harmonic oscillator 
bath. This implies that the density matrix is Gaussian so 
that environmental information is contained in the sec- 



ond moments (q 2 ) and (p 2 ) [EQ], 



(qW) 



exp 
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2(q 2 ) 
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(6) 

Expectation values of higher powers of H s are non-trivial 
because q and p do not commute. The purity of the 




FIG. 2: The probability to measure a harmonic oscillator in 
the ground and first three excited states as a function of x and 
y (see text). The line traces out the behavior of the ohmic 
bath as a function of the coupling in the under-damped range. 



density matrix Eq. JBJl is 



Trp 2 = J dqdq'{q\p\q'){q'\p\q) = 



k/2 



x/WW 2 



(7) 



The uncertainty relation, */ (q 2 ) (p 2 ) > ft/2, guarantees 
that Trp 2 < 1, with the inequality becoming sharp if 
the oscillator is isolated from the environment. As the 
environment causes greater deviation from the Planck 
scale limit, the state loses purity. 

The generating function Z may be calculated conve- 
niently by tracing in the position basis and inserting a 
complete set of position states between the operators, 



Z(x)= / dqdq'{q\p\q'){q'\e-^\q). 



(8) 



The first object in Eq. JSJ is the density matrix in posi- 
tion representation, given by Eq. ©. The second object 
may be interpreted as the uncoupled position-space prop- 
agator of the harmonic oscillator from position q to q' in 
time —ihx- We find 



2E 



sinh e\ 



2 A (coshex — 1) 



1 + cosh e\ 



where e — fko. 2E 



(9) 

\q 2 ) + (p 2 }/m and A = 
(q 2 )(p 2 )/h 2 . E is the average energy of the oscillator, 
while A > 1 is a measure of satisfaction of the uncer- 
tainty principle. Eq. has a pleasing limit for the free 
particle u> — > 0, 



Z(x)free = {l + x(p 2 )/m} 



(10) 



which is just the generating function for Wick contrac- 
tions, (p 2n ) = (2n- l)!!((p 2 ))". Thus, in Eq. 10, the 
inverse square root generates the right combinatorial fac- 
tors under differentiation, and the nontrivial x depen- 
dence accounts for the commutation relations between q 
and p. The first few harmonic oscillator energy cumu- 
lants may now be straightforwardly found via Eq. (JJl , 



«# s 2 » 
((# s 4 » 



(l/2)[-(e 2 /2) + 4E 2 -2e 2 A] , 
-(E/2)[~16E 2 + e 2 (l + 12A)] 
48£ 4 - 4e 2 £ 2 (l + 12,4) 
e 4 [(l/8) + 2A + 6A 2 }. 



(11) 
(12) 

(13) 



After inserting the mean square values for an ohmic bath 
(see the discussion above eqs. H15ll6| l). Eq. (|llf> is iden- 
tical to the main result of Ref. [171 ]. 

Alternatively, we now consider the diagonal matrix 
elements p nn . An analytical expression for the den- 
sity matrix in the energy basis may be found by using 
the wavefunctions of the harmonic oscillator, ip n {q) cx 
e -7 i / 2 H n (^q) where 7 = yjmuj/h and H n (x) is the n th 
Hermite polynomial. In the energy basis, the density ma- 
trix is given by p nm = / dqdq'ip* n {q){q\p\q')^ m {q'). The 
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position space integrals may be done using two differ- 
ent copies of the generating function for the Hermite 
polynomials. The diagonal elements may be found by 
equating equal powers of the generating variables. We 
first define the dimensionless variables x = 2 r y 2 (q 2 ), y — 
2{p 2 ) /{^ 2 h 2 ), and D = 1 + x + y + xy. x and y are related 
to the major and minor axes of an uncertainty ellipse. 
The isolated harmonic oscillator (in it's ground state) 
obeys two important properties: minimum uncertainty 
(in position and momentum) and equipartition of energy 
between average kinetic and potential energies. The in- 
fluence of the environment causes deviations from these 
ideal behaviors which may be accounted for by introduc- 
ing two new parameters, a = (y — x)/D, b — (xy — 1)/D 
with — 1 < a < 1 and < b < 1. The deviation from 
equipartition of energy is measured by a, while the devi- 
ation from the ideal uncertainty relation is measured by 
b. We find 



reveals entanglement properties of the many-body quan- 
tum mechanical ground state. Consequently, repeated 
experiments on simple quantum systems give information 
about the nature of the environment, the strength of the 
coupling and entanglement. The larger the energy fluc- 
tuations, the greater the entanglement. There are sev- 
eral possibilities for experimental implementations. We 
have mentioned measurement of persistent current 0, @ 
as well as projecting on the system's energy eigenstates. 
Another measurement possibility is a zero temperature 
activation-like process |l8| where the dominant mecha- 
nism is not tunneling, but the same quantum effects of 
the environment which we have discussed here. 

This work was supported by the Swiss National Science 
Foundation. 




m+2p • 



(14) 



m— p— 



All terms are positive (because a appears only squared), 
less than one, and obey Trp = 1. The summation over 
m and p define a polynomial P n of order n in a and b. 
The probability for the lone oscillator to be measured in 
an excited state clearly decays rapidly with level num- 
ber. These probabilities also reveal environmental infor- 
mation. For example, P\ — 2b and is thus only sensitive 
to the area of the state, while Pi — a 2 + 2b 2 depends 
on both the uncertainty and energy asymmetry. Addi- 
tionally, if we expand the first density matrix eigenvalue 
0, El w ith respect to small deviations of x and y, we 
recover pu in agreement with our general argument. 

Although x and y have been treated as independent 
variables, the kind of environment the system is coupled 
to replaces these variables with two functions of the cou- 
pling constant. For example, with the ohmic bath fljlH^ ] 
(in the under-damped limit), the variables are 



x{a) 



2 a 
1 arctan , 

7T 



y(a) = (1 - 2a 2 )x(a) + — In — 



. (15) 
(16) 



where a is the coupling to the environment in units of 
the oscillator frequency and uj c is a high frequency cutoff. 
This bath information is shown in Fig. J2J) with uj c = lOw. 
The trajectory of the line over the surface shows how the 
probabilities evolve as the coupling a is increased from 
to 1. Other kinds of environments would trace out 
different contours on the probability surface. 

In conclusion, we have shown that projective measure- 
ments of the system Hamiltonian at zero temperature 
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